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Abstract: Resource-based competition between microorganisms species in continuous culture has 
been studied extensively both experimentally and theoretically, mostly for bacteria through Monod 
and and Contois "constant yield" models, or for phytoplankton through the Droop "variable yield" 
models. For homogeneous populations of N bacterial species (Monod) or N phytoplanktonic species 
(Droop), with one limiting substrate and under constant controls, the theoretical studies [TJ[21[3] 
indicated that competitive exclusion occurs: only one species wins the competition and displaces 
all the others. The winning species expected from theory is the one with the lowest "substrate 
subsistence concentration" s* , such that its corresponding equilibrium growth rate is equal to the 
dilution rate D. This theoretical result was validated experimentally with phytoplankton [4] and 
bacteria [5], and observed in a lake with microalgae [BJ. On the contrary for attached bacterial 
species described by a Contois model, theory [7] predicts coexistence between several species. In 
this paper we present a generalization of these results by studying a competition between three 
different types of microorganisms : free bacteria (represented by a generalized Monod mode), 
attached bacteria (represented by a Contois model) and free phytoplankton (represented by a Droop 
model). We prove that the outcome of the competition is a coexistence between several attached 
bacterial species with a free species of bacteria or phytoplankton, all the other free species being 
washed out. This demonstration is based mainly on the study of the substrate concentration's 
evolution caused by competition; it converges towards the lowest subsistence concentration s*. 
leading to three different types of competition outcome: 1. only the free bacteria/ phytoplankton 
best competitor excludes all other species; 2. only some attached bacterial species coexist in the 
chemostat; 3. A coexistence between the best free species, with one or several attached species. 

Key-words: competition, competitive exclusion, droop, variable yield model, monod, ratio- 
dependent, biomass-dependent, microorganism, microalgae, phytoplankton 



Competition entre phytoplankton et bacteries: 
exclusion et coexistence 



Resume : La competition pour la ressource entre micro-organismes dans des 
cultures en continu a ete largement etudiee experimentalement et theoriquemcnt , 
surtout entre bacteries modelisees par des taux de croissance de type Monod ou 
Contois, ou pour le phytoplancton a travers le modele de Droop. Pour les pop- 
ulations homogenes composees de N especes bacteriennes (Monod) ou N especes 
phytoplanctoniques (Droop), avec un substrat limit ant et des commandes main- 
tenues constantes, les etudes theoriques (TJ [3] ont indique que l'exclusion com- 
petitive se produit: une seul espece remporte la competition et elimine toutes les 
autres. L 'espece dont la theorie predit la victoire est celle avec la concentration s* 
de substrat permettant sa survie la plus basse. Ce resultat theorique a ete validee 
experimentalement pour le phytoplancton [4j et les bacteries [5], et a ete observe 
dans un lac avec des microalgues [BJ. Par contre, pour les especes bacteriennes 
decrites par un modele Contois, la theorie predit la coexistence entre plusieurs es- 
peces (7]. Dans cet article nous presentons une generalisation de ces resultats en 
etudiant une competition entre les trois differents types de micro-organismes: des 
bacteries libres (representees par un modele de Monod generalise), des bacteries 
fixees (representees par un modele de Contois) et du phytoplancton (represente par 
un modele de Droop). Nous prouvons que le resultat de la competition est une 
coexistence entre plusieurs especes bacteriennes fixees avec une espece de bacteries 
libres ou de phytoplancton, toutes les autres especes libres etant lessivees. Notre 
demonstration est basee principalement sur l'etude de 1'evolution du substrat causee 
par la competition; elle converge vers la plus faible concentration de subsistance s*, 
ce qui conduit a trois types differents des resultats de la competition: 1. seule la 
meilleur bacterie libre ou le meilleur phytoplancton exclut toutes les autres especes; 
2. seules quelques especes bacteriennes fixees coexistent dans le chemostat, 3. une 
coexistence entre la meilleure espece libre, avec une ou plusieurs especes fixees. 

Mots-cles : competition, exclusion competitive, droop, monod, ratio-dependence, 
microorganisme, microalgues, phytoplancton 
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1 Introduction 

1.1 Growth of phytoplankton 

Phytoplankton is composed of microscopic plants at the basis of the aquatic trophic 
chains. Phytoplankton means a broad variety of species (more than 200.000) using 
solar light to grow through photosynthesis. Phytoplankton plays a crucial role in 
nature since it is the point from which energy and carbon enter in the food web. But 
it may also be used in the future for food or biofuel production, since several phy- 
toplankton species turn out to have very interesting properties in terms of protein 
[5J[9] or lipid [10 content. In addition to light, phytoplankton requires nutrients for 
its growth. The "paradox" of phytoplantkon species coexistence was introduced by 
Hutchinson [TT|: "The problem that is presented by the phytoplankton is essentially 
how it is possible for a number of species to coexist in a relatively isotropic or un- 
structered environment all competing for the same sort of materials". In this paper 
we consider this question from a theoretical viewpoint "what are the mechanisms 
leading to competitive exclusion or coexistence, and to what competition outcome 
do they lead?". But so far most of the competitions studies have assumed that only 
phytoplankton species were engaged in the competition. However, it is clear that 
such species also have to compete with the bacteria for nutrients. 

In this paper, we study the competition between phytoplankton and bacteria. 
Phytoplankton can be accurately represented by a Droop model |12( [T3J [H] which 
accounts for their ability to store nutrients and to uncouple uptake and growth. 
Bacteria are represented by simpler models. They can be of two different types, 
described either by a Monod type model if they live in suspension or by a contois 
model if they are attached to a support. 

In this paper we first recall the main results available for competition of microbial 
species of the same class. Then we consider the problem of 3 class competition. After 
some mathematical preliminaries we state and demonstrate our main Theorem. A 
discussion concludes our paper and highlights the ecological consequences of our 
result. 

1.2 The Competitive Exclusion Principle (CEP) 

"Complete competitors cannot coexist" 

This is the formulation chosen by Hardin [JS] to describe the Competitive Exclusion 
Principle (CEP). According to him, this ambiguous wording "is least likely to hide 
the fact that we still do not comprehend the exact limits of the principle". But still, 
a more precise formulation is given: if several non-interbreeding populations "do 
the same thing" (they occupy the same ecological niche in Elton's sense |T5]) and 
if they occupy the same geographic territory, then ultimately the most competitive 
species will completely displace the others, which will become extinct. 

Darwin was already expressing this principle when he spoke about natural se- 
lection ([17] p. 71 and 102). Scriven described and analyzed his work in these words: 
"Darwin's success lay in his empirical, case by case, demonstration that recogniz- 
able fitness was very often associated with survival. [...] Its great commitment and 
its profound illumination are to be found in its application to the lengthening past, 
not the distant future: in the tasks of explanation, not in those of prediction" |18j . 
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Since the work of Darwin, men have tried to apprehend the limits of the principle 
in different context and by different means. In the next sections we present how 
mathematical models have shown their appropriateness for predicting the outcome 
of competition, in the case of chemostat-controlled microcosms. 

1.3 The chemostat, a tool for studying the CEP 

"Microbial systems are good models for understanding ecological processes at all 
scales of biological organization, from genes to ecosystems" [IS]. The chemostat is 
a device which enables to grow microorganisms under highly controlled conditions. 
It consists of an open reactor crossed by a flow of water, where nourishing nutrients 
are provided by the input flow, whereas both nutrients and microorganisms are 
evacuated by the output flow. To keep a constant volume in the vessel, these two 
flows are kept equal. In this paper we consider that the following conditions are 
imposed in the chemostat: the medium is well mixed (homogeneous); only one 
substrate is limiting for all the species, whose only (indirect) interaction is the 
substrate uptake; the environmental conditions (temperature, pH, light, ...) are 
kept constant, and so are the dilution rate D, corresponding to the input/output 
flow of water, and the input substrate concentration Sj„. Figure [l] represents such 
a chemostat. 



Sin 




Figure 1: A chemostat, which enables to grow microorganisms under highly con- 
trolled conditions. The input/output flow of water is D, and the input substrate 
concentration is Sj„ 



The chemostat has been used to study the CEP since the beginning of the XX 
century |20| . and its experimental use has often been coupled with mathematical 
models [2]. 

1.4 Bacterial and phytoplanktonic models, and previous the- 
oretical results on single class competition 

1.4.1 Free bacteria growth (generalized Monod model) 

To predict the growth of bacteria in suspension within a chemostat, Monod devel- 
oped a model [21], where the growth rates of the biomasses xi (i G {1, • • • , N x } for 
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a competition between N x species) depend on the extracellular substrate concentra- 
tion s. In the classical Monod model the growth rates cti(s) are Michaelis-Menten 
functions 

<*(«) = TT^r 

where a™ are the maximum growth rates in substrate replete conditions, and Kf 
are the half saturation constants. In this paper we consider a generalized Monod 
model to represent growth of free bacteria, by using the wider class of functions 
verifying Hypothesis [T] 

Hypothesis 1 M-model: 

a,-(s) are C , increasing and bounded functions such that aj(0) = 0. 
We note a™ the supremum of the growth rate: 

supers) = a" 1 > 

The free bacteria dynamics write 

±i = (ati(s) - D)xi 

with s, x % € R+ for i e {1, • • • , N x ] and D € M+. 1 J 

In this model the substrate uptake is proportional to the biomass growth for each 
bacterial species, so that the total substrate uptake per time unit will be a i( s )^ 

1.4.2 Phytoplankton model (generalized Droop model) 

Phytoplankton is able to uncouple substrate uptake of nutrients from the growth 
associated to photosynthesis |13| . This capacity to store nutrients can provide a 
competitive advantage for the cells that can develop in situations where substrate 
and light (necessary for phytoplankton growth) are rarely available concomitantly. 
This behaviour results in varying intracellular nutrient quota: it is the proportion 
of assimilated substrate per unit of biomass Zk\ it can be expressed for instance 
in mg[substrate]/mg[biomass]. Droop [32] developed a model where these internal 
quotas are represented by new dynamic variables qk (denoted "cell quota"). The 
substrate uptake rates Pk(s) are assumed to depend on the extracellular substrate 
while the biomass growth rates Jk(<lk) depend on the corresponding cell quota. 

In the classical Droop model the functions have specified forms. Uptake rates 
are Michaelis-Menten functions ^ of the substrate concentration: 

Pk{s) = JTlq pT (2) 

and the growth rates are Droop functions (J3j) of the cell quotas: 

= { (3) 
1 if q k < Q° k 

with pj™ and jk the maximal uptake and growth rates; represent the half satura- 
tion constants, and Q° k the minimal cell quota. In this paper we consider the wider 
class of Q-models (Quota models) verifying Hypothesis [2j so that it can encompass, 
among others, the classical Droop formulation [T2] as well as the Caperon-Meyer 
model 1221. 
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Hypothesis 2 Q-model: 

• pk(s) are C , increasing and bounded functions such that p k (0) = 

• Ikilk) are C > increasing and bounded functions for q k > > 0. When 
Qk < Ql> IkiQk) = 0. 

It directly ensues from Hypothesis [2] that fk(qk) = 7k(qk)qk are increasing functions 
(for qk > Ql) which are onto R+, so that the inverse functions f^ 1 are defined on 

We denote p™ and ~/ k the supremal uptake and growth rates: 

sup s > Pfc(s) = pf > 
su Pg k >Q ( l lk(qk) = Ik > 

The phytoplankton dynamics write 
<lk = Pk(s) - fk(qk) 

ik = (ik(qk) - D)z k (4) 
with s, q k , z k e M. + for k € {1, • • • , N z } and D e K+. 

The substrate uptake per time unit is Ylk=i Pk(s) z k 

This model has been experimentally shown to be better suited for phytoplank- 
ton dynamic modelling than the Monod model ([2]) that implicitly supposes that 
the intracellular quota is simply proportional to the substrate concentration in the 
medium and which must be definitely limited to bacterial modelling. The stability 
of the Q-model has been extensively studied in the mono-specific case ( (2"51 |2"H 123] ). 

1.5 Previous demonstrations of the CEP for M- and Q-models 

The advantage of Monod and Droop models is that their relative simplicity allows a 
mathematical analysis. The analyses of the M-model with N x bacterial competing 
species [T], and of the Droop model with 2 phytoplankton species [2] and then 
recently with N z phytoplankton species [3] led to a confirmation of the CEP in the 
chemostat, and to a prediction on "who wins the competition", or "what criterion 
should a species optimize to be a good competitor". In both cases, we have 

Theorem 1.1 // environmental conditions are kept constant and the competition 
is not controlled (D and Si n remain constant) in a chemostat, then the species with 
lowest "substrate subsistence concentration" sf* (or s k *), such that its corresponding 
equilibrium growth rate is equal to the dilution rate D, is the most competitive and 
displaces all the others. 

A striking point about this result is that it permits to make predictions on the result 
of a competition, only by a priori knowledge of the species substrate subsistence 
concentrations sf* (or s z k *). This latter can be determined in monospecific-culture 
chemostat, so that the competition outcome can be determined before competition 
really occurs. Several experimental validations where carried out with phytoplank- 
ton [4] and bacteria [5]. This theoretical behaviour was also confirmed in a lake [6], 
where the species with lowest phosphate or silicate subsistence concentrations won 
the competition for phosphate or silicate limitations. 
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1.5.1 Attached bacteria model (generalized Contois model) 

In case where bacteria are not free in the medium but there is a spatial heterogeneity 
{e.g. they grow attached on a support, such as floes in the culture medium), a ratio- 
dependent model is more adapted to describe bacterial growth. Contois model f261j 
represents such dynamics by using more complex growth functions where the growth 
rates depends on the ratio of the substrate concentration over biomass concentration 
.'/.. (j G {I,-" ,N y }): 



fr(«.yj) = WxTrf 



In this paper we consider the wider class of "C-model" (Contois model), wich is 
more general than a ratio dependent model. It verifies the following hypotheses: 

Hypothesis 3 C-model : 

f3j(s,yj) are C 1 functions on K + X K + \ {(0,0)}, increasing and bounded functions 
°f s (f or Uj > 0), and decreasing functions of yj (for Sj > 0) such that Vt/j G 
M+, 3j(0, yj) = and Vs G K+, lim pAs, yA = 

We also need to add the following technical hyposthesis, which is verified by the 
classical Contois function: 



Hypothesis 4 



d 



We notice that the Contois growth function is undefined in (0,0), and that Hypoth- 
esis^has been built so that this property can (but does not have to) be retained by 
the generalized j3 function. All other properties imposed by Hypotheses^ and^ are 
satisfied by the original Contois growth-rate. 

We denote Pj l {yj) the supremal growth rates for biomass concentration yj: 

supp j (s,yj)=pj l (y j ) 

s>0 

so that the C-species dynamics write 

Vj = Wj{s,Vj) - D) yj , , 

with s, yj G R + for j G {1, ■ • ■ , N y } and D G E+. K ' 

In this model, like in the M-model, the substrate/biomass intracellular quotas bj 
are supposed to be constant for each species, so that the substrate uptake rates are 
proportional to the growth rates with a factor 1/bj. 

1.5.2 Coexistence result for competition between C-species 

Competition between several C-species was studied fTjj and led to a coexistence at 
equilibrium with the substrate at a level s v * depending on the input substrate concen- 
tration Si n and the dilution rate D. The species share the available substrate. To be 
more precise we must define the ' 'sq- compliance" concept: SQ-compliant species are 
the species able to have a growth rate equal to the dilution rate D with a substrate 
concentration so- The results offTjj show that all the "s y * -compliant" species coexist 
in the reactor at equilibirum, and all the others are washed out, as they cannot grow 
fast enough with substrate concentration s y * . 
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Definition 1.2 A species Xi,yj or Zk is sq- compliant if it is able to reach a growth 
rate equal to the dilution rate D with a substrate concentration sq. 

1.5.3 Competition and coexistence - towards a new paradigm 

Following these results an interrogation arises: 

■C What would be the result of a competition between "competitive " free bacterial 
and microalgal species, and "coexistive" attached bacterial species? Competitive 

exclusion? Coexistence? 3> 

The aim of this paper is to provide an answer to this question, and to give insight 
into the mechanisms forcing the outcome of such a competition. This answer leads 
to a broader view and understanding of competitive exclusion and coexistence mech- 
anisms, following the words of Hardin f!5p "To assert the truth of the competitive 
exclusion principle is not to say that nature is and always must be, everywhere, "red 
in tooth and claw. " Rather, it is to point out that every instance of apparent co- 
existence must be accounted for. Out of the study of all such instances will come 
a fuller knowledge of the many prosthetic devices of coexistence, each with its own 
costs and its own benefits. " 

1.6 A generalized model for competition between several 
phytoplankton and bacteria species growing according 
to different kinetic models 

The generalized model for competition between all bacteria and phytoplankton species 
is an aggregation of these models, which alltogether give the following substrate 
dynamics, subject to substrate input, output, and uptake rates: 

N y N z 

s = D(s in - s) -^ai(s)— -^2fij(s,yj)^ -J2p k (s)z k (6) 
i=i ai j=i °i fc=i 

The parameters related to the nutrient flow are the dilution rate D > and the 
input substrate concentration s$ n > 0, which are both assumed to be constant. 

To simplify notations we can remark that this system can be normalized with 
ai = bj = 1, when considering the change of variables Xi = ^ and yj — |^ (note 
that all the hypotheses are still satisfied). We obtain system where variables Xi 
and yj are now expressed in substrate units. 

s = D(s in - s) -22ai(s)xi ^ \, (•-•// .i//. ~^2pk(s)z k 

i = l j = l k=l 

< ±i = (a,(s) - D)xi 
Vj = (ft ( s ,%) - D )Vj 

Zk = (ik(qk) - D)z k 

. q\k = pk(s) - fk(qk) 

with fkiik) = ik(qk)qk 

ands,q k <E R + , s in ,D G M+.and Xi(0), yj(0), Zk{0) € for i e {1, 
j G {!,••• ,N y },ke {!,••• ,N Z } 



■■■,N X }, 
(7) 
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Note that the results obtained in this paper apply also on the simple M- only, Q- 
only, and C'-only competition models, or on a model with two of these three kind of 
species. 

1.7 Other coexistence mechanisms, and competition control 

This introduction wouldn't be complete without a short review of what has been done 
concerning other coexistive models, or the control of competition. 

Following the question arised by Hutchinson JTTjj concerning the "paradox of the 
phytoplankton ", a large amount of work has been done to explore the mechanisms 
that enable coexistence, mainly for models derived from the Monod model. It has 
been shown to occur in multi-resource models ]27[ \28\j . in case of non instantaneous 
growth f29\/ . in some turbidity operating conditions ]3(J^ . a crowding effect or 
variable yield ]3°2?$ (not in the Droop sense). ]33^ and \3$ also presented several 
mechanisms which can mitigate the competition between microorganisms and pro- 
mote coexistence. 

In other papers ()33$. f36\/ and 1311), controls were proposed to "struggle against 
the struggle for existence" (that is, to enable the coexistence of complete competi- 
tors). These controls indicate how to vary the environmental conditions in order to 
prevent the CEP from holding : some time varying or state- depending environmental 
conditions can enable coexistence. f38\j propose a theoretical way of driving compe- 
tition, that is, of choosing environmental conditions for which the competitiveness 
criterion changes. 

2 Mathematical preliminaries 
2.1 The variables are all bounded 

Throughout this paper we study the evolution of one solution of system with 
initial condition (s(0), cci(0), . . . ,2^(0), 

I/i(0),...,yjv y (0),g 1 (0),...,g JV ,(0),«i(0),... 1 */v.(0)) where x t (0) > 0,^(0) > 0, 
Zfc(0) > 0. In this section we study the boundedness of the variables. First, the 
variables all stay in M. + , as their dynamics are non negative when the variable is 
null. 

Then we know that the biomasses remain positive: 
Lemma 2.1 

V* G {1, • • • , N x }, xi{0) > <S> Vt, Xi(t) > 
Vj 6 {1, • • • , N y }, yj (0) >Q<^Vt, Vj (t) > 
Vfce {l,-- - 7 N z },z k (0) > 0^Vt,z k (t) > 

Proof: Because of the lower bounds on the dynamics (±i > —Dxi for the free 
bacteria for example), the biomasses are lower bounded by exponentials decreasing 
at a rate D: 

\/t,Xi(t) > Xi{0)e- Dt > 

□ 
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Then, to upperbound the variables we define 

JV, N y N z 

i=l j=l k=l 

the total concentration of intra and extracellular substrate in the chemostat. The 
computation of its dynamics gives 

M = D(s in - M) (8) 

so that M converges exponentially towards s in . This linear convergence implies the 
upper boundedness of M: 

yt > 0, M{t) < M m = max(M(0), s in ) 

Then s, Xi,y.j and q k z k are also upper bounded: 

Vi > 0, s(t) < M(t) < M m 

Vie{l,--- ,JV x },Vt>0, Xl {t)<M{t)<M m 

Vj e {1, • • • , N y }, Vt > 0, y 3 (t) < M(t) < M m w 

Vfce{i,--- ,N z }yt>o, qk {t)z k {t) < M(t) < M m 

We are now interested in the boundedness of the Q-model's cell quotas q k and 
biomasses z k 

Lemma 2.2 Vfc, the q k variables are upper bounded by max(f^ (ph(M m )) , gfc (0) ) 

Proof: For any q k > f k ~ 1 (p k (M m )) there is an upper bound on q k : 

Qk = Pk(s) - f k (q k ) < p k {s) - Pk {M m ) < 
so that s < M m implies that q k cannot increase if it is higher than fj7 (pk(M m )). □ 



Lemma 2.3 Vfc, the z k variables are upper bounded by 

( M m \ 
zr = max ^r— ,z k (0)) (10) 

with the convention that 7 / 7 1 (Z?) = +oo if ~f k < D 

Proof: As q k z k is upper bounded by M m , there is an upper bound on z k 

z k = (7fe(9fe) - D)z k < ^7 fe (^—^J ~ z k 
so that z k cannot increase if it is larger than ^r^y- □ 



Lemma 2-4 After a finite time to there exists a lower bound § > for s. 
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Proof: With hypothesis and as the biomasses are upper bounded, we see that s 
can be lower bounded 

s > D(s in - s) - g - £ &(«, y?)y? - £ - 0(a) 

i=i j=i fc=i 

where cj> is a decreasing function of s, with 0(0) = Dsi„ and 0(s,- n ) < 0. -By continu- 
ity of the 4> function, there exists a positive value s < Si n such that 4>{s) — Dsi n /2. 
The region where s > s is therefore positively invariant. Also s is increasing for any 
value lower than s with s > Dsi n /2 so that s(t) reaches s after some finite time to- □ 



Remark 1 This lemma eliminates any problem that could have arisen from, the 
problem of definition of /3j(s,yj) in (0,0). After the finite time to, no solution can 
approach this critical value anymore. 



Lemma 2.5 There exists a finite time t\ > such that for any time t > t\, 
quit) e (Qg, QT) mth QT = fkHpT)- 



Proof: If qk(t) > QT > then we have 

Qk < Pk(s) - fk{Qk) < Pk(M m ) ~ pT < o 

for all gfc G [QT i <7fc(0)], so that qk(t) < Q™ in finite time t± and for any t>t±. 
U Qk(t) < Ql with t > to (defined in Lemma 2.4-), then we have that 

qk = Pk(s) > p k (s) > 

for all qj, G [qk(to), Ql], so that qk(t) > * n finite time t\ and for any t>t\. 



□ 



This lemma is biologically relevant since minimum and maximum cell quotas are 
indeed known characteritics of phytoplankton species. For the rest of this paper we 
will consider that all the q^ are in the {Q° k ,QT) intervals. 

Remark 2 In the classical case of Michaelis-Menten uptake rates and Droop 
growth rates |#|) we have: 

m 

QT = Ql + ^ 

ik 



2.2 Prom a "substrate" point of view... (How substrate con- 
centration influences the system) 

Since model is of dimension 1 + N x + N y + 2N Z , it is hard to handle directly. 
In this section we introduce functions which clarify how the qk and yj dynamics 
are influenced by s. This will enable us to focus on the substrate concentration 
evolution, and thus reduce the dimension in which the system needs to be analyzed. 
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2.2.1 Internal cell quotas qu are driven by the substrate concentration 

s 

It is convenient to introduce the functions 

Q k (s) = f,; x (p k (s)) (ii) 

and 

S z k (q k )=Qk\<lk) (12) 

With Hypothesis^it is easy to check that Qk is defined, continuous, increasing from 
(0, +oo) to {Q Q k , Q™), so that S k is also well defined, continuous and increasing from 
{QIjQT) t° (Oj+oo). The qk equation can then be written 

q k = fk(Qk(s)) - fk(q k ) (13) 



or 

q k = Pk(s) - Pk{S z k (q k )) (14) 

Since fk(qk) and pk{s) are increasing functions, we see how the dymanics of q k is 
influenced by the sign of Qk(s) — qk (or s — S k (qk)): 

sign{q k ) = sign(Q k (s) - q k ) = sign(s - S k (q k )) (15) 

For a given constant substrate concentration s, the equilibrium value of qk is Qk(s). 
Conversely, s must be equal to S k (qk) for q k to be at equilibirum. 

Function Qk realizes a mapping from the substrate axis to the cell quota axis. 
Functions S k realizes a mapping from the cell quota axis to the substrate axis. An 
illustration of the cell quotas behaviour is presented in Figure^ 

91 Ql{s) 

e — * • *— e 



Figure 2: Two equivalent statements: "<j% goes towards Q k (s)" and "S%.(qk) goes 
towards s" (see the sign Property (15)). The latter permits a one dimensional view 
of the s and qk dynamics, on the substrate axis. 



2.2.2 How the biomasses yj are driven by the substrate concentration s 

For the C-species, it is also convenient to introduce functions Yj(s): 

if f3j(s,0) > D, then Yj(s) is defined by /3j(s,Yj(s)) = D , . 

if Pj{s, 0) < D, then Yj(s) = ' 

and the inverse Sj(yj) functions: 

w ^ n / ^ 3s ° s± ft( s °'%) > D -> then Sj(yj) is de fi ned b v P( s j(yj)>yj) = D 

V% >U '\ else, S]( yj )=+oo 
S?(0)=w£ yj>0 S]( yj ) 

(17) 
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The values of s such that Yj(s) — correspond to values where the substrate is 
too low for yj to survive (yj is not s-compliant at these values). The values of yj 
such that Sjfyj) = +00 correspond to levels of biomass yj that cannot be sustained 
independently of the substrate level. 

With Hypothesis [5] it is easy to check that Yj is defined, continuous, increasing 

from (S V AQ), +00) to I 0, sup Yj(s) I , so that is also well defined, continuous and 
\ s >a J 

increasing from I 0, sup Yj(s) I to (S*(0), +00) . 

The yj equation can then be written 

y j = (p j (s,y j )-p j (s,Y j (s)))y j (18) 

or 

Vj = (fis(*,Vi)-0 3 iS*ivj),Vi))vi (is) 



Thus with yj positivity (see Lemma 2.1) we see how the dymanics of yj are influenced 
by the sign ofYj(s) - y (or s - Sjljjj)): 

sign(yj) = sign(Yj(s) - y 3 ) = sign(s - (20) 

For a given constant substrate concentration s, the equilibrium value of yj is Yj(s). 
Conversely, s must be equal to Sj(yj) for yj to be at equilibirum. 

Function Yj realizes a mapping from the substrate axis to the cell quota axis. 
Functions Sj realizes a mapping from the cell quota axis to the substrate axis. An 
illustration of the biomasses behaviour is presented in Figure [5[ 

2/1 Yxis) 

Y-2{s) y-2 
* • -^3 



Figure 3: Two other equivalent statements: "yj goes towards Yj(s)" and "S v Ayj) 



goes towards s" (see the sign Property (20 1). The latter permits a one dimensional 



view of the s and yj dynamics, on the substrate axis. 



Finally, with Figures ^| and [5] we obtain a one dimensional view of the s, qk 
and yj dynamics on the substrate axis. The demonstration presented in this paper 
ensues mainly from this one dimensional view of the system. 

2.3 The convergence of s is related to the convergence of qk 
and yj 

Lemma 2.6 In system the five following properties are equivalent for any sq > 
mi nj -(5}(0)): 
i) lim t ^ +oc s(t) = s 
ii) yi,\im t ^ +00 q k (t) = Q k (s ) 
Hi) 3i,\im t ^ +00 q k (t) = Q k (s ) 

iv) Vj, lim^ +00 yj(t) = Yj(s ) 

v) 3j, lim t ^ +oc Uj(t) = Yj(s ) > 
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When lim t _j. +oc s(t) = sq < min J (S'J(0)) ; all the qk(t) converge to Qk(so) an d the 
y 3 (t) to Yj(s o ) = 0. 



Proof: In the case sq > min J (5j'(0)) we successively demonstrate five implica- 
tions. 

i => ii and i => iv: straightforward with the attraction (13) of q^ by Qk{s), and 
the attraction (18) of yj by Yj(s). Note that yj(0) cannot be null (Lemma 2.1). 

ii => Hi and iv => v: trivial implications. 

Hi => i (and v => i): we equivalently demonstrate that the simultaneous conver- 
gence of qk (resp. yj) and non convergence of s lead to a contradiction. 



\s(t) - s 




Figure 4: Visual explanation of the demonstration of Lemma |2.6| s is repeatedly 
escaping a ^-interval around sq (•)■ Because \s\ is upper bounded by B, then s is 
out of the ry/2-interval during non negligible time intervals (dashed lines represent 
|s| = B). qk (resp. yj) is repeatedly attracted away from Qk(so) (resp. Yj(so)) by 
Qk{s) (resp. Yj(s)) (arrows) 



If s does not converge towards Sq, it is repeatedly out of a [sq — r\, Sq + n] interval, 
denoted r/-interval: 

3r) > 0,Vi > 0,3t s > t, \s(t s ) - s | > r] 
In Figure^ t s time instants are represented by 



We can then use the upper-bounds and (10) on s, xi, yj and z\, to show the 
boundedness of the s dynamics 



D(s m -M m )~J2 a 



mm 
i J 'i 



i=l 



N y N z 

g /3f (0)yf - £ P ^4 n < s < Ds in (21) 
i=i 



fc=i 



so that 



\s\<B 



with B = max | Ds in , -D(s in - M m ) + a™x\ n + ^ /3f (0)y? + ^ - 

\ i=l j=l k=l 

Then, every time s is out of the rj-interval, it must also have been out of the r//2- 
interval during a time interval of minimal duration A(i]) — (For a visual 
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explanation see the dashed lines of Figure representing the increase caused by 
\s\=B). 

If for some t s we have s{t s ) > s + r\, we then have that s(t B ) > So + T]/2 during 
the whole time-interval [t s — A(rj),t s ]. We can thus lower bound the dynamics of qk 
(resp. yj) during that time-interval: 

qk = fk(Qk(s)) - fk(qk) 

> fk(Qk(s + r,/2))-Mq k ) 

and 

Vj > (ft Oo + 17/2, yj) - D)yj 

= (ft(s + n/2, Vj ) - ft(s + rj/2, Yjiso + f]/2))) Vj 

Now the convergence of q k to Q k (so) (resp. yj to Yj(so)) is defined as 

Ve > 0,3t« > 0,Vt > t\ \q k {t) - Q fc (so)] < e{resp. \y,(t) - Y 3 (s )\ < e) (22) 

since we can pick e such that e < Q(sq + rj/A) — Q(so) (resp. e < min(Yy(so + 
?7/4) — Yj(so), Yj(so) — Yj(so — ry/4)) ), we then have, for t > t q , that qk(t) < Q(sq + 
rj/A) (resp. Yj(sQ — rj/A) < yj{t) < Yj(sq + rj/A)). Taking our t s larger than the 
corresponding t q + A(q), we then have for all time t € [t s — A(q),t a ] 

q k > fk(Qk(s + T?/2)) - (fk(Qk(so) + e)) 

> fk(Qk(so + - fk(Qk(s + »?/4)) = C«(r]) > 

and 

Vj > (ft ( s o + t}/2, Yj (s ) + e) - ft (s + ij/2, Y - (s + ?//2)))% 

> (ft (s + ?7/2, ^-(so + ri/4)) - ft (s + ?7/2, Yj(8 + v/2)))Vj 

> (ft (s + 77/2, 5^ (s + r?/4)) - ft (s + t?/2, Y}( flo + tj/2)))^ (s - 77/4) 
= Cv(rf) > 

C q (q) > since is an increasing function of qk, and C v (j]) > since ft is 
a decreasing function of yj 
We then define 

C(r?) =min(C«(77),C ,2 '(77)) 

// we t/ien choose e such that e < MlLh^M ; ^/j en |g fc (7j s ) — gj.(t s — A(ry)) | > 2e 
(resp. \yj(t s ) — yj(t s — A(r)))\ > 2e), and we see that the increase of qk (resp. yj) 
makes it eventually get out of the e-interval around Qk{so) (resp. Yj(so)). This is 
a contradiction, so that implication iii => i (resp. v => i) holds. 

Alternatively, if for some t s we have s(t s ) < s — r), then we can upper bound 
the dynamics of qk (resp. yj) during the [t s — A(r)),t s ] time-interval: 

Qk < fk(Qk(so~ V/ 2 )) - fk(Qk) 
and 

Vj < (ft («o ~ v/ 2 , Vj) ~ D) Vj 

= (ft (s - n/2, Vj ) - ft (s - 77/2, Y,-(s - r}/2))) Vj 

and the same arguments hold, with 

q k < fk(Qk(s - ??/2)) - f k (Qk(s ) - t?/4) = C( V ) < 
and 

yj < (ft (s - 77/2, Yj(s - 77/4)) - ft (s - 77/2, Yj(s - 77/2)))^ (s + 77/4) 
= C( V ) < 

Inria 



Competition between phytoplankton and bacteria 



17 



and finally e < C^YM*!) yjfacfo causes the contradiction. □ 



2.4 The equilibria correspond to the substrate subsistence 
concentrations 

In this section we present the equilibria of the generalized competition model 
The first equilibrium of this model corresponds to the extinction of all the microor- 
ganisms species: 

E = (s in , 0, ... ,0 , 0, ... ,0 , 0, ... ,0 , Qi(s in ), . . . ,Qjv(sm)) 

This equilibrium is globally attractive if the input substrate concentration Sj„ is not 
high enough for the species' growth to compensate their withdrawal of the chemo- 
stat by the output flow D, that is if \fi, cti(si n ) < D and Vj, /3j(s,- n , 0) < D and 
Vfc, 7fe(Qfe(sm)) < D (proof of this result is easy and we omit it; for getting a clear 
idea of the demonstration, see JTjj and l^j for the Monly and Q-only cases). We 
suppose that we are not in this situation through the following hypothesis: 

Hypothesis 5 We assume that one of the following condition is satisfied: 

• 3i, ai(s ln ) > D 

. 3j,p j (s in ,0)>D 

• 3k,j k (Q k (s in )) > D 

This guarantees that, at least for one of the families of species, there exists some 
index i, j, k and some associated unique sf*,s^ , s^* < s in (denoted "subsistence 
concentration") such that 



7*(0*(*D) = D 



P j {s y j *,Y j (sf))=D with sf + Y 3 {sp 



Note that in the C-model, there exists an infinity of s € [5j(0),Sj„) verifying 
(3j(s,Yj(s)) = D. The value s y * is then the substrate concentration required for 
having species j remaining alone in the chemostat at equilibrium. It has to satisfy 
s|* + Yj(sJ*) = Si n because of that imposes M = s + yj = Si n at equilibrium. 
We number these species such that 



U < s^ = sf* < 8%" < . . . < s~* < s in 

o < sm < s y 2 (o) < ■ ■ ■ < sih < Sin 

< s z * = sf < sf < ... < s z n * < s in 

with V(i,i, k) e ({1, ...,n x }, {1, ...,n y }, {1, ...,n z }), sf* ? S](0) ? s%* 



(23) 



where n x , n y and n z are the number or free bacteria, attached bacteria and phyto- 
plankton species having a subsistence concentration smaller than Si n for the given 
D; all other species cannot be positive at equilibrium. Hypothesis [3| implies that at 
least one of n x , n y and n z is non-zero. We denote s x * and s z * the lowest M- and 
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Q- substrate subsistence concentrations. We also denote s y * the substrate concen- 
tration that there would be at equilibrium if there were only attached species in the 
chemostat (see J?fl/); since it needs to satisfy |^p, it requires 

n y 
3=1 

Though the sum ofYj(s v *) spans all the relevant indices, some species might have 
Yj(s v *) = because they have s y * < S y (0) < s^ n . If some n x , n y or n z is zero, we 
set the corresponding s x * , s y * or s z * to Si n because none of the species from their 
family can survive at a substrate concentration lower than Si n , which is the higher 
admissible concentration. 

In the previous competitive exclusion studies JIJ O \7jl these quantities were of 
primer importance, as they directed the result of competition. Here we show that 
the competition outcome is strongly linked to 

s* = min(s x *,s v *,s z *) 

which is the lowest of all subsistence concentrations. Hypothesis [3| implies that 
s <C s^ n . 

We do not consider the case where two subsistence concentrations are equal, 
because we suppose that the biological parameters of each species are different. In 
his broad historical review about competitive exclusion ]13$ Hardin wrote: "no two 
things or processes in a real world are precisely equal. In a competition for substrate, 
no difference in growth rate or subsistence quota can be so slight as to be neglected". 



Hypothesis 6 V(i, j, k) G ({1, n x }, {1, nj, {1, ...,n z }), sf* ^ S»(0) + s%* 

The subsistence concentrations and Yj(s) functions are presented in Figure [5[ 
In this figure, we see that no free species model can coexist at equilibrium because 
s cannot simultaneously be equal to sf* and s^* . On the contrary, attached species 
verifying Hypothesis |5| can support different s value at equilibrium (between 5j(0) 
and Si n ), so that there exist equilibria where one M- or D-model species coexist with 
one or several attached species (see Figure^for a graphical explanation). On those 
equilibria, only the attached species verifying 

5V(0) < sf (resp. 5j(0) < s^) (24) 

can coexist as they can be at equilibrium at the subsistence concentration of the M- 
model (resp. D-model) species, by having a biomass equal to Yj(sf*) (resp. Yj(s^*)). 

For these considerations, we can enunciate the following proposition which does 
not need to be proved: 

Lemma 2. 7 For a given sq substrate concentration, we have 

• Xi is so-compliant if sf* = So; 

• yj is so-compliant is Sj(0) < so; 

• Zk is so-compliant if s" = sq 
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Figure 5: Subsistence concentrations of the M-model (sf*) and Q-model (si*) 
species, and equilibrium biomass Y, (s) of the attached species, which enable these 
species to have a growth rate equal to the dilution rate D 7 and thus to be at equilib- 
rium. We see that M- and Q-model species cannot coexist at equilibrium because 
they have only one fixed subsistence concentration, and s cannot be simultaneously 
equal to several of these concentrations. On the contrary, attached species can co- 
exist with others at equilibrium because they can have a growth equal to D for any 
s £ [Sj(0),Si n ), by adjusting their biomass concentration to Yj(s) (see definition 
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It ensues that, for the corresponding C-species we have 

Y j (s )>0 

which means that they can be at positive equilibrium under dilution rate D and 
substrate concentration Sq. Thus, the s* -compliant species are: 

• only the s* -compliant C-species, if s* = s y * 

• x\ and all the s* -compliant C-species, if s* = s x * 

• Z\ and all the s* -compliant C-species, if s* = s z * 

We now present all these equilibria and their stability in M-, C- and Q-only substrate 
competitions. 

2.5 M-only equilibria 

Ef = (sf , r; i) , 0,...,0 , 0,...,0 , Q x («f* ),..., QnM*)) 

with X* — Si n — sf* 

each of these M-only equilibria corresponds to the winning of competition by free 
bacteria species i; such an equilibrium only exists for i £ {0, • • • ,n x } (all other 
species cannot survive at a substrate level lower than Si n fot the given D). In 
a competition between several free bacteria, eqilibrium Ef (with lowest substrate 
subsistence concentration sf* ) is asymptotically globally stable, while all the others 
are unstable JJ]/. 

2.6 Q-only equilibria 

E z k = {sf , 0,.„,0, 0,...,0, 0,...,4,...,0, Qi(sI*),...,QnM*)) 

Similarly to M-only equilibria, each of these phytoplankton only equilibria correspond 
to the winning of competition by phytoplankton species k; such an equilibrium only 
exists for k £ {0, • • • , n z }. In a competition between several phytoplankton species, 
equilibrium E\ (with lowest substrate subsistence concentration s\* ) is asymptoti- 
cally globally stable, while all the others are unstable J^j. 



2.7 C-only equilibria 

We denominate G a subset o/{l, ••• ,n y } representing any C-species coexistence. 
For example if we want to speak about species 1, 5 and 7 coexistence, then we use 
G = {1,5,7}. We then define E V G the equilibrium where these species coexist. It is 
composed by 

• s V q such that s V q + X)jeG ^i( s G*) = Sin because of 

. VjeG, % = y,-( s «*) 

• for any other j , yj = 

• Vie {!,••■ ,N x },Xi = 
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• VfcG {I,-- - ,N z },z k =0 

. Vfce{l,-- - ,N s }),q k =Q k (s%) 

there exist many E G equilibria, corresponding to all the possible G subset. The 
globally asymptotically stable equilibrium of a competition with only attached species 
is given by the choice G — {1, • • • ,n y } f^j/. Note that some of the G species can 
have a null biomass on these equilibria, as Yj(s G ) might be null for some j € G. 
Therefore E G and E Go with G\ ^ G2 are not necessarily different. 

We must here introduce a technical hypothesis which will be useful later to prove 
hyperbolicity of the equilibria. 

Hypothesis 7 For all G and all j : Sj(0) 7^ s v G 



2.8 Coexistence equilibria 

As previously said in this section, there also exist equilibria where one of the free 
species coexist with several sf*- or sj?* -compliant attached bacterial species. For a 
coexistence with free bacteria species we denote them E^ X q V ^* . They are composed 
of: 



• for any other j , yj = 

• V/ ^ i, xi = 

. Vfce {l,... ,N z },z k = 

. Vfce{l,-- - ,N z },q k =Q k (sf) 

• Xi = Si n — sf* — ~^2j£G^j( s f*) (this value will be denoted xf ) 

Similarly, for a coexistence with phytoplankton species we denote them E^ G ^* 
They are composed of: 

• s = 4* 

• for any other j , yj = 

. V*e{l,--- ,N x }, Xi = 
»Vl^k, Zl =0 

. Vie{l,---,JV,},gj = Q I («f) 

• z k = — — Sfc qfff*') Y ^ Sk - (this value will be denoted z^ ) 
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To our knowledge, these equilibria have never been studied until now. 

Note that some of these equilibria might be reduntant with M- or Q-only equilib- 
ria, if all the C-species represented by G are not s x * - or s|* -compliant. Note also 
that all those equilibria do not necessarily exist in the non-negative orthant. In- 
deed, xf and zj? can be negative, depending on Si n and on the substrate subsistence 
concentrations. These equilibria with negative components will not be studied any 
further since we only consider initial conditions in the positive orthant, which is 
invariant. In the sequel, we will denote E an equilibrium of which belongs to an 
unspecified class. 

We will now show that if s — s* at equilibrium, there exists a positive equilibrium 
containing all s* -compliant species. 

Lemma 2.8 

• If s* — s x * then E^^ _* n j is in the positive orthant. 

• If s* — s v * then E^* n j is in the positive orthant. 

• If s* — s z * then E^^* m n j is in the positive orthant. 

Proof: 

• If s* = s x * , then all yj — Yj(s x *) > at equilibrium and s v * + Y^j=i Yj(s y *) = 
Si n implies that s x * + X)j=i Yj(s x *) < s in since s x * < s v * and Yj(s) is non- 
decreasing. It directly follows that x\ = s in — s x * — Yj{s x *) > 0. 

• If s* = s y * , then all Xj and z% are zero at equilibrium and all yj — Yj(s y *) > 

• If a* = s z \ then all y. } = Y^s 2 *) > at equilibrium and s v * + E"=i Y j(s y *) = 
Si n implies that s z * + Y^,j=i Yj{s z *) < Si n since s z * < s v * . It directly follows 

s, n -s"-Y,™ z Yj(s z *) 

that z\ = Qk(s**) > °- 

□ 

We call E* the equilibrium with all s* -compliant species remaining in the chemo- 
stat, while all the others are excluded. Depending on the species subsistence concen- 
trations, E* can be one of the previously presented equilibria: 

• if s* — s y * then E* — E y * n , : only the s* -compliant C-species remain in 
the chemostat. 

• if s* — s x * then E* = E^f^* n , : the best free bacteria species (lowest s x * ) 
remains in the chemostat with all the sf* -compliant C-species. 

• if s* = s z * then E* = E^ z &* n j : the best phytoplankton species (lowest si* ) 
remains in the chemostat with all the sf,* -compliant C-species. 

In the next section, we present an important global stability result for this equilib- 
rium. 
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3 Statement and demonstration of the Main Theo- 
rem: competitive exclusion or coexistence in the 
generalized competition model 

This theorem states that all the s* -compliant species (those who can be at equilibrium 
with substrate subistence concentration s* , which is the lowest of all s x *.s y *,s z *) 
coexist in the chemostat at equilibrium, while all the others are excluded. 

Main Theorem 1 In the generalized competition model if Hypotheses 

hold, then all the solutions of the system, having CEj(O), yj[0), £fc(0) > for all 
s* -compliant species, converge asymptotically towards equilibrium E* . 

Structure of the proof: In a first step we reduce system to the mass balance 
surface. Then we present a non decreasing lower bound L(t) for s(t), and use it to 
demonstrate that s converges towards s* . Finally only the s* -compliant species have 
a large enough substrate concentration to remain in the chemostat, so that all other 
M-, Q-, and C-model species are washed out. The final step consists in showing that 
the convergence result that we showed on the mass-balance surface can be extended 
to the whole non-negative orthant. 

Remark 3 It is not restrictive to consider Xi(0),yj(0), £fe(0) > for the solutions 
of the system since species with null initial condition can be ignored, so that we can 
then consider a smaller dimensional system. 



3.1 Step 1: we consider the system on the mass balance 
surface and in the region where E (Q%, Q™) for all 
k e{l,...,N z } 



Lemma 2.5 indicates that q k reaches {Q kl Q™) in finite time, and in |<5|) we showed 
that the total concentration of intra and extracellular substrate in the chemostat M 
converges to Si n . 

We denote "S ", the generalized competition model on the mass balance sur- 
face defined by 

JVx N v N z 

m = s + Xj + yj + q k z k = s in (25) 
1=1 j=i fe=i 

For the remainder of the demonstration we will study system E, and we will later 
show that its asymptotic convergence towards an equilibrium has the same behaviour 
as the initial model While studying system £, we will however retain all the 

states of the original system and the expressions of the equilibria; E is then defined 
by the addition of the invariant constraint (25). 



3.2 Step 2: we propose a non decreasing lower bound L(t) 
for s 

The main obstacle for the demonstration of the Main Theorem was the possibil- 
ity that s would repeatedly be lower than s* and repeatedly be higher than = 
max(s",s^*), which would generate an oscillating behaviour. In order to eliminate 
this possibility we build a non decreasing lower bound for s, which converges towards 
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s*. We now present such a lower bound, which will be used to show that s converges 
to s* in the next sections. 

Lemma 3.1 In system E 

L(t) = min (min(S z k (q k {t))), wm(S» (%(*)), a*,a(t) 
\ fc j 

is a non decreasing lower bound for s 

Proof: We know that the right derivative of L is the derivative of one of the 
function which realizes the minimum. In four cases we show that this right derivative 
is non negative. 

• Case 1: If S k {q k {t)) realizes the minimum then its derivative is non negative, 
because St(q k ) goes towards s (see (15)). 

• Case 2: If Sj(yj(t)) realizes the minimum then its derivative is non negative, 
because Sj(s) goes towards s (see (20)). 

• Case 3: If s(t) realizes the minimum then we examine its dynamics s for 
system E (i.e. on the mass balance equilibrium manifold). We replace Si n by 
s + Ei x i + Ej Vj + Efe 1kZ k : 

s = ^2(D - oti(s))Xi + ^2(D - Pj(s,yj))yj + ^{Dq k - p k (s)) z k 

i j k 

which is equivalent to, from the definition ofQk(s) : 

S = ^2( D - a i( s )) x i + ^2( D - ft'( s ) Vj))Vj + ^2 ( Dqk ~ 7k(Qk(s))Qk(s)) Z k 

i j k 

Then 

— for all i, s < s* gives us aj(s) < D, so that the first sum is non negative; 

— for all j, s < S v 3 {y )) gives us f3j(s,yj) < Pj{S^{yj),yj) = D), so that 
the second sum is non negative; 

— for all k, s < S k (qk) gives us Qk(s) < q k , and s < s* gives us Jk(Qk( s )) ^ 
7fc(Qfe(Sfe*)) = D so that the third sum is also non negative. 

Finally we obtain 

s > 

• Case 4 : If s * realizes the minimum, we know that its right derivative is null 
and thus non negative. 

□ 



3.3 Step 3: we demonstrate that s converges towards s* 
Lemma 3.2 In system E 

lim s(t) = s* 

t— v+oo 
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Proof: We first show, by contradiction, that the substrate concentration s(t) can- 
not converge towards any constant value other than s* . Suppose the reverse hy- 
pothesis, i.e. lini t _j. +00 s(t) = s ^ s* . Through Lemma 2.6 we then have that 
lim t ^ +00 q k (t) = Qfe(s) and lim t _j. +00 y 3 (t) = Y 3 (s). 
Ifs< s*, 

• cti(s) < D for all i so that all Xi go to 

• 7fe(Qfc(s)) < D for all k implies that all z k go to 



So that we have a contradiction with mass balance equilibrium {25), as the total 
substrate (in the medium + in the biomasses) at equilibrium s 

+ E"ZiY 3 (s) will be 

lower than s y * + Xw=i Yj(s v *) = Si n . 

If s > s* we must consider three cases: 

• if s* = sf* then Qfj(s) > D implies that x\ diverges to +oo, which is in 
contradiction with the boundedness shown in 

• if s* — sf* then 7i(Qi(s)) > D implies that Z\ diverges to +oo, which is in 
contradiction with the boundedness shown in (10). 



• if s* — s v * then we have a contradiction with mass balance equilibrium (25), 
because s + J^i=i ^j'(^) w iH ^ e higher than s y * + Y]j—i Yj(s v *) < Sj n . 

Hence the impossibility of convergence of s towards any s other that s* is proven. 
We now demonstrate the lemma by contradiction. We assume that 

s does not converge towards s* 

which, from the previous remark means that s does not converge to any constant 
value. 

Remark 4 As s does not converge towards s* , we know that the q k do not converge 
towards Qk(s*) (see Lemma 2.6) 

We consider two cases, which both lead to a contradiction, on the basis of a reason- 
ing which is close to the demonstration developed to prove Lemma \2.6\ 

• Case a: L attains s* in finite time 

In Appendix^A^we show that a contradiction occurs. 

Idea : s cannot stay higher than s* without converging to s* , because this would 
cause x\ or z\ to diverge, or s + X^i^K s ) ^° be always higher than Si n without 
converging to Si n . 

• Case b: L never attains s* 
See Appendix\B\ 

Idea : If s did not converge to s* , the non decrease of L and its attraction by s 
would cause it to reach s* . 

In both cases we found a contradiction, so that the proof of Lemma \3.2\ is com- 
plete. 

□ 
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3.4 Step 4: all the s*-compliant species remain in the chemo- 
stat, while the others are excluded 

In this section we show that, as s converges towards s* in model S, all the free 
species with substrate subsistence concentration higher than s* are washed out of 
the chemostat because their growth a.i(s) or "fk(Qk) cannot stay high enough to com- 
pensate the output dilution rate D. Finally, all the s* -compliant species able to be 
at equilibrium with a substrate concentration s* remain in the chemostat. 



Lemma 3.3 In system £ all the solutions with positive initial conditions for the 
s* -compliant species converge to E* . 



Proof: For all the xi and Zk species such that on(s*) < D and 7fe(Qfe(s*)) < D, it 
is straigthforward that the convergence of s to s* will cause their biomass to converge 
to 0. If s* — s y * , then this is true for all the free species. 

For all the s* -compliant C-species, we have from Lemma \2.6\ that their biomass 
will tend to Yj(s*), which is positive for the s* -compliant species and null for all 
the others. 

Finally, if s* = s x * or s* = s z * , then we have through the mass balance equi- 



librium (25) that the free species whose subsistence concentration is s* will have its 

— N y 
biomass converge to Si n — s* — z2j=i Yj( s *) : a ^ substrate which is not present 

in the medium or in the attached biomasses is used by the best M- or Q-competitor. 

□ 



3.5 Step 5: convergence of the solutions for model £ implies 
convergence for model ([7]) 

In order to extend the convergence result to the full model and thus prove our Main 
Theorem, we apply a classical theorem for asymptotically autonomous system ]3<A l^j. 

Lemma 3.4 All solutions of system with positive initial conditions for the s* - 



compliant species converge to E* defined in section 2.7 



Remark 5 While, up to here, we simply had considered £ as the same system as 
|?]) ; in the same dimension, except that it was restricted to (25), we will now equiv- 



alent^ explicitely include (25) into system to obtain £ in the form of a system 
that has one dimension less than by omitting the s coordinate. Since both repre- 
sentations of £ are equivalent, the previously proven stability results are still valid 
in the new representation, with the exception that convergence takes place towards 
equilibria directly derived from these presented in sections \2.5§2.8\ by omitting the s 
coordinate. These new equilibria are differentiated from the original ones by adding 
a ~ , so that an arbitrary equilibrium is denoted E. 
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Proof: System S can be written as follows : 



a, 



E 

1=1 



Xl 



y.i 



E x < 

1=1 



Ny 

E f" 

m— 1 

£ Vn 

m—1 



E 

r=l 

E 

r=l 



q r z r , yj 



D au 



D Vj 



(26) 



h = hk(qk) - D)z k 



qk = pk 
for i e {l, • 



N r 



E 

m— 1 



, a;/ - y ,y m - 7 . q r z r 
,N X }, j €{!,-■■ ,N v },k€{l,- 



E< 

r=l 



- fk(qk) 
,N Z } 



where the s state has been removed compared to |?|). In order to 
we should add the equation 

M = D{s m - M) 



recover 



model 



which we interconnect with (26) by replacing every Sj„ in (26) with M. It is this 
interconnection that we will now study. 

In the first part of the proof, we will show that every solution of converges 
to an equilibrium E. We will then show by induction that all the solutions that do 
not converge to E* have an initial condition with some Xi = 0, yj = or z k = for 
some s* -compliant species. Thus, all the solutions with Xi 7^ 0,yj 0, z k 7^ for 
the s* -compliant species converge to E* . 

For that, we will use Theorem F.I from J^j. We will therefore first compute the 
stable manifolds of all equilibria o/S: 

• The stable manifold of E* is of dimension N x + N y + 2N Z . It is constituted of 
all the initial conditions which verify ^(0), Vj{0), Zfe(0) > for s* -compliant 
species and Xj(0), y 7 (0 ), Zfc(O) > for all other species, as well as q k > for 
all k (see Lemma 3.3). 

• The stable manifold of Eq is of dimension N x — n x + N y — n y + 2N Z — n z . It 
is constitued of all the initial conditions which verify X\(0) = ... = a;„ x (0) = 
0, yi(0) = ... = y ny (0) = and Zl (0) = . . . = z„ z (0) = 0. The only 
species that can be present at the initial condition are those that cannot survive 
for the given D and Sj„. Indeed, if any a^(0) > for i < n x (or similar 
yj(0) > or Zfe(0) > 0), one can apply Lemma 3.3 to the reduced order 
system containing these species to show that convergence does not take place 
towards Eq. Conversely, any initial condition with £Ei(0) = . .. = x nx (0) = 0, 
2/i(0) = ... = 2M y (0) = and Zi(0) = ... = ^(0) = generates a solution 
that goes to E since for the other species we have: 

— Xi < (ai(si n ) — D)xi, with cti(si„) — D < for all i > n x because of the 
definition of n x presented in (23); 

— yj < (/3j(Sj n , 2/j) — D)yj, with /3j(sj n ,0) — D < for all j > n y , because 
of the definition of n y ; 

— 4 < (lk(Qk(sin)) - D)z k , with j k (Qk(sin)) - D) < for all k > n z 
because of the definition of n z . 
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The di mens ion of the stable manifold of any other E can be computed from 
To an equilibrium E corresponds a substrate value s (> s* by 



3.3 



Lemma 

definition of s* ). Lemma 3.3 indicates that solutions o/E converge towards an 
equilibrium corresponding to s, if there is no smaller subsistance concentration 
corresponding to a species present in the system (for free species) and if all 
M-, Q- and C-species that are s-compliant are present in the corresponding 
equilibrium. The stable manifold of E must therefore be constrained to initial 
conditions that verify ajj(0) = 0, Dj(0) = and Zfc(0) = for all species that 
are s-compliant for some s < s and that are not positive in E. Having set all 
these values to zero, it is indeed clear that s is the s* as defined in Lemma 



3.3 of the reduced order system (without the aforementionned Xi,yj and Zk 



coordinates). All solutions defined in Lemma 3.3 of this system then converge 
to E, which justifies our definition of the stable manifold of E. Its dimension 
is N x + N y + 2N Z — n E s , where n E ~ is the number of s-compliant species (for 
some s < s) that are not present in E. 



Through Lemma \3.3\ we have in fact shown that all solutions of E in the non- 
negative orthant converge to an equilibrium. Indeed, for a given initial condition, 
either it belongs to the stable manifold of E$ or, eliminating from the system all 
species that are null at the initial time necessarily sets it in a form where Lemma 



3.3 can be applied (which shows convergence to an equilibrium). 

The dimension of the stable manifold of any equilibrium E will therefore be the 
one of E plus 1. The hypotheses of Theorem F.l from [21 are indeed all verified: 

• The whole system is bounded (see section 2.1) 



The equilibria of system E are hyperbolic (see Appendix C.2\C6 ). 



• There are no cycles of equilibria in system E. Indeed, if we analyze the poten- 
tial transition between two equilibria, both equilibria must belong to the same 
face, so that convergence takes place to the one corresponding to the smallest 
value of s. A potential sequence of equilibria would then be characterized by a 
decreasing value of s at each equilibrium, which prevents it from cycling. 

We can then conclude from this theorem that all solutions of tend to an equilib- 
rium. We are then left with checking to what equilibrium they tend. 

Before continuing this proof, we need to detail n E ~. In the case of E — E* and 
s = s* , we have s = (by definition, all s* -compliant species are present in E* 
and there is no other species that is compliant for smaller values of s). Otherwise, 
we necessarily have n E ~ > 0. Indeed, we know that s > s* , so that all species 

present in E* are compliant for some s < s; as such, in order to have n E - = 0, E 

would need to at least contain all species that are present in E* . In such a case no 
S and Q species can be present in E* ( otherwise, it could not be present in E also 
for a different value of s). Defining J the set of C-species that are present in E* 



and writing (25) for E* then yields 



M = s * + J2y^ 
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Equality (25) should also be valid in s > s* so that 



N y N z 

Sin = S + ^ X, + ^ V] + ^2 °-k Z k > S + ^ Y 3 ( S ) > S in 

i=l j=l k=l j&J 

where we have the last inequality (which leads to a contradiction) because Yj(s) is 
an increasing function. We can then conclude that, for all E ^ E* , - > 0, and 
at least one s* -compliant species species must be null. 

In order to check to what equilibrium solutions of |?|) tend, we use an induction 
argument, by supposing that our Main Theorem has been proven up to N — 1 species, 
which we use for the proof for N species. Along with the fact that the stability result 
is trivial for 1 species (classical Monod model, [2}, classical Droop model, \2J$ and 
generalized Contois model, J3j), this will conclude our proof. 

Let us consider a system of N species with equilibrium E* as defined earlier. 
This equilibrium contains positive species (which are s* -compliant) and null species 
(which are not s* -compliant). 

Imposing, for one of the not s* -compliant species, Xi = (or yj — or Z], = 0) 
for the initial condition, sets us in the framework where we have N—l species present 
in the system. Also, since this species did not belong to the positive ones in E* , its 
absence does not change anything into which equilibrium is the one corresponding 
to the smallest subsistance concentration, which remains E*. We can then apply 
the induction hypothesis, which indicates that all such initial conditions initiate 
solutions that converge to E* (as long as the s* -compliant species have positive 
initial condition). 

Studying now the equilibrium Eq, we know from the beginning of the proof that 
its stable manifold is of dimension N x — n x +N y — n y + 2N z — n z + l. As was done for 
S, it is directly apparent that any initial condition with a:i(0) = . . . = x ni .(0) = 0, 
t/i(0) = . . . = y n (0) = and z±(0) = . . . = z rix (0) = generates a solution that has 
all species exponentially go to zero. Finally, the analysis of the s equation shows 
that it has the form s = D(si n — s) — F(t) with F(t) exponentially going to zero so 
that s goes to Si n and all such solutions go to Eq. 

We can now consider all the other equilibria. Let an equilibrium E corresponding 
to a substrate concentration s (> s* by definition) . As we have seen in our analysis 
of S, the stable manifold of the corresponding E is of dimension N x + N y + 2N Z — 
so that the stable manifold of E is of dimension N x +N y + 2N Z —n^ -+ 1. Let 
us set ourselfes in the situation where all s species are set to zero at the initial 
time and all others are positive. We can then consider the system with only the 
remaining N x + N y + 2N Z — ~ positive species and the substrate. We have seen 
that, in this case, all solutions of the corresponding reduced order E go to E which 
means that s is the "s*" defined in Lemma \3.J\ for the reduced order system. Since 
the reduced order system contains less than N species because - > 0, we conclude 
that all solutions of the full system that have zero initial condition for all - 
species and positive values for all N x + N y + 2N Z — ~ others converge to E. We 
have then exhibited an invariant manifold of dimension N x + N y + 2N Z — - + 1 
for which all solutions go to E; this corresponds to the predicted dimension of the 
stable manifoldof E. No solution with some of the n^- species positive (among 
which there is at least on s* -compliant species) at the initial time can then converge 
to E. 
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This completes the proof of our Main Theorem since all solutions go to an equilib- 
rium and we have exhibited the stable manifold of all equilibria other than E* . These 
manifolds cannot go into the region where X{,yj or Zk > for all s * —compliant 
species because at least one of them is in the corresponding ns s -set. All initial con- 
ditions in the region where Xi,yj or Zk > for all s * —compliant species therefore 
generate solutions that go to E* . □ 



4 Discussion 

4.1 How D and s in both determine competition outcome 

In M- and Q-only competitions, the outcome of competition is mainly determined 
by D, which fixes the sf* and sf* M- and Q-substrate subsistence concentrations; 
the role of Si„ is to allow the best competitor (already determined by the value of 
D) to settle the reactor, or to cause it to be washed out with all the others. On 
the contrary in C-only competition, both controls have important roles: D fixes the 
Yj(s) functions, while Si n determines the equilibrium, where s y * Yj(s v *) — Sj„. 

With a low enough Sj n , only few C-species will settle the chemostat (s y * being low 
in this case, there will be few s y * -compliant species, with non-null Yj(s v *)), whereas 
a high enough Sj„ can enable all C-species to coexist. 

Finally, in a mixed competition the dilution rate D fixes all the M- and Q- 
substrate subsistence concentrations sf* and sf* , as well as the Yj(s) functions, while 
the input substrate concentration Si n selects the species remaining in the reactor, by 
limiting the available nutrients, and thus the biomasses present in the reactor at 
equilibrium. Figure [#| gives an example between three competitors. 

On this figure the s x * and s z * values and the Y(s) function are fixed by D. 
Here s\* is lower than sf* , so that the free bacteria species will be outcompeted and 
washed out. Then the value of Sj„ determines wether 

1. no species remain at equilibrium 

2. only the attached bacterial species remains at equilibrium, as there is not 
enough input substrate to feed both attached bacteria and phytoplankton species: 
because s y * + Y(s y *) = Si n and Si n < s z * + Y(s z *), we know that s y * < s z * , 
so that s* = s v * , and only the C-species remains in the reactor. 

3. both the attached bacteria and phytoplankton species remain in the chemostat: 
here s m > s z * + Y(s z *) and s ln = s v * + Y(s v *) give s y * > s z * , so that 
s* = s z * and the phytoplankton species remains in the reactor, coexisting with 
the s z * -compliant C-species. 

In this last case D has fixed the s z * substrate equilibrium value and the Y(s) func- 
tion, and at equilibrium the total substrate in the chemostat, equal to Si n , will be 
composed of 

• the substrate in the medium s z * (which is fixed by D and does not depend on 

Sin) ' f 

• the attached bacterial species internal substrate Y{s z *) (which is also fixed by 
D only); 
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Figure 6: Mixed 3 class competition outcome depends both on the dilution rate D 
and input substrate concentration Si n . The solid lines represent the influence of D, 
which fixes the subsistence concentrations of one M-model (s x *) and Q-model (s z *) 
species, and the equilibrium biomass Y(s) of one C-species. As the free bacteria 
species has a too high subsistence concentration s x * > s z * it will be outcompeted 
and excluded. The three numerated zones represent the influence of Si n . Zone 1 
(si 7l < S v (0)) : no species remain at equilibrium. Zone 2 (S v (0) < Si n < s z * + 
Y(s z *)) : only the C-species remains at equilibrium. Zone 3 (sj„ > s z * + Y(s zi ')) : 
the attached bacteria and phytoplankton species coexist. 



• the phytoplankton species internal substrate Q(s z *))z* = Si n — Y(s z *) — s z * , 
which depends on Si n . 

By going from left to right in Figure [#| starting with Sj n = 0, it is possible to 
imagine the input substrate concentration increase, thus enabeling more and more 
substrate s — Si n at equilibrium (zone 1). Then in zone 2 the C-species is present 
at equilibrium, and as Si n increases, more and more biomass Y{s) is present at 
equilibrium. Finally Y(s z *) is the maximal biomass for which the attached species 
needs less substrate at equilibrium than the phytoplankton species to have a growth 
rate equal to D. After that it has to coexist with the phytoplankton species: when s in 
increases higher than s z * + Y(s z *) it enables more and more Q-biomass z* , while 
keeping substrate concentration s — s z * and C-biomass y = Y(s z *). 



4.2 Originiality of the demonstration 

The demonstration explains how the state variables evolve, and its originality for 
the study of uniquely phytoplankton (or bacteria) species can be summed up in three 
points. 

First, we chose to study the substrate evolution instead of ignoring it after the 
classical mass balance equilibrium transformation s — Si n — x i~^2j Vj~^2k IkZk- 
Then, the definition of the Sj and Sf, functions enabled to gather most information 
on the substrate axis: instead of having separate information on 1 + N x + N y + 2N Z 
axes we obtained a one dimensional view on these dynamics (Figure \B and \3y, 
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where all the Sj(yj) and S^(qk) go towards s. We have thus turned a complex 
1 + N x + N y + 2N Z dimensional problem into a simpler one: "how do s and the 
Sj(yj) and S^(qk) behave on the substrate axis, and what are the consequences for 
the biomasses?". 

Finally the definition of the non decreasing lower bound L(t) (section 3.2) and its 
convergence towards s* (section 3.3) were the last steps for this demonstration to 
emerge. 

Free species pure competitions (with one class of species among Monod or Droop) 
for substrate lead to the "survival of the fittest ", the fittest being the species with 
lowest substrate requirement s* . On the contrary, Contois-only competition lead to a 
coexistence equilibrium, because biomass dependence gives attached bacterial species 
the capability to remain at equilibrium for different substrate concentrations in the 
range [5^(0), Sj n ) (see Figure^. Monod and Droop species are mutually exclusive, 
which leads to the pessimization principle of adaptative dynamics 140/ ■' "mutation 
and natural selection lead to a deterioration of the environmental condition, a Ver- 
lenderung. We end up with the worst of all possible environment. " On the contrary 
attached species are coexistence- compliant thanks to biomass dependence, which nu- 
ances the pessimization principle: "some species could live in worse environments 
(s = mhij(Sj(0)) being the worse one) but if there is enough substrate for other 
species, they can coexist." (see Figure^and discussion) 



5 Conclusion 

In this paper a demonstration was given for the outcome of competition between 
phytoplankton and bacteria. Three scenarios are possible, depending both on the 
dilution rate D and input substrate concentration Si n (see discussion for precisions): 

• only the best free competitor remains in the chemostat; 

• only some attached bacterial species coexist at equilibrium; 

• a new equilibrium (never studied before) is attained, where the best free com- 
petitor coexists with all the s* -compliant attached bacterial species. 

Since the introduction of the concept of evolution, with its link to competitive 
exclusion f!5^ and the "paradox of phytoplankton " modelling has tried to ap- 
prehend competition, and to predict or control it. Our contribution in this frame- 
work was to extend the results proven in the N-species Monod model, N-species 
Droop model and N-species Contois model, where the outcome of competition was 
predicted and explained with mathematical arguments, accompanied by ecological in- 
terpretations. 

An important conclusion in this type of competition is that attached bacteria are 
likely to be present in a pure culture of microalgae. This may have very important 
consequences on the ecological point of view, since such natural coexistence between 
a phytoplanktonic species and attached bacteria may have lead to co-evolution, where 
the best association between phytoplankton and bacteria have been progressively se- 
lected. 
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A Step 3 - Case a: L attains s* in finite time 

In this case 

• if s* = sf* we consider Figure [?] where L attains s* after a finite time t L : 

V< > t L ,L{t) = s* 



Qi(*!) + 7/2 



Figure 7: Visual explanation of the demonstration of Lemma |3.2| - Case 1: L attains 
s* in finite time t L (Q-model). i) qi is repeatedly higher than Qi(s*) + 8 (•). ii) 
Because qi is upper bounded by p™, so that qi is higher than Qi(s*) + 9/2 during 
non negligible time intervals (dashed lines represent q\ = p™). Thus Z\ diverges, 
which is a contradiction. 




Substep 3a.l: after a finite time larger than t L , q± is repeatedly 
higher than Qi(s*) + 8. 

Since mhii(S%.(qk)) > L, we know that 

\ft>t L , qi {t) >Qi(s*) 

As s does not converge to s*, we also know from Lemma \2.6\ that qi does not 
converge towards Qi(s*): 

39 >0,\/t> 0,3t q > t, \ qi {t q ) - Qi(s*)\ > 8 

Those two facts imply that the repeated exits of qi(t) from the 8-interval around 
Qi(s*) take place above Qi(s*) for any t q > t L , so that, in that case, we have 
q\(t q ) > <2i(s*) +8. In Figure^ such t q time instants are represented by •. 

Substep 3a.2: q\ is higher than Qi(s*) + 8/2 during non negligible 
time intervals. 

Since the qi-dynamics are upper bounded with 

91 < P? 
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we know that every time q\ is higher than Qi(s*) + 0, it has been higher than 
Qi(s*) + 0/2 during a time interval of minimal duration A(8) = On 
Figure [?j q\ — p™ is represented by the dashed lines. 



Substep 3a.3: then z\ diverges, which is impossible 

From time t L on, we have that q\ > Qi(s*) =>■ 7i(<7i) > D, so that z±(t) 
is non decreasing. During each of the time interval where q\ is higher than 
Qi(s*) + 0/2, the increase of z\ is lower bounded by 

*i = 7i(0i(**) + 0/2) ~D = C(8) > 
so that every t q time we have 

Zl {t q ) ~ Zl {t q - A{6)) > C{6)A{6) 
As such increases occurs repeatedly, and as Z\ is non decreasing, Z\ diverges. 



This is a contradiction because Z\ is upper bounded (see {10)). 



• if s* = sf* then the non convergence of s to s* , and the fact that s > s* will 
cause s to be non negligibly "away" from s* , so that x\ will diverge, causing 
a contradiction with This is exactly the same demonstration as above (in 
the case s* = sf*) without needing the q^ study. 

• if s* — s v * then s + 5Zfci^j( s ) ™^ always be higher than s in = s v * + 



Ei=i y j( sV *) without converging to Si n , which is in contradiction with 125V 



B Step 3 - Case b: L never attains s* 

In this case L(t) converges towards a value L £ (0, s*], because it is non decreasing 
and bounded in [0, s*], so that 

Ve>0,3t L (e) > 0, Vi > t L {e), \L{t) - L\ <e 

We consider the neighborhood of L in Figure\£\ 



Substep 3b. 1: after a finite time, s is repeatedly higher than L + X. 



Since, from the beginning of the proof of Lemma \3.S\ we know that s does not 
converge to any constant value, hence not to L, 

3X > 0, V<> 0, 3t s > t, \s(t s ) - L\ > A 

Since L is increasing and converges to L, it reaches L — A in finite time t L (X). 
After this finite time, s is higher than L + X on every t s time instants, which are 
represented by • in Figure[S\ 
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Figure 8: Visual explanation of the demonstration of Lemma |3.2| - Case 2: L 
never attains s*. i) s is repeatedly higher than L + A (•). it) s is upper bounded 
by Dsi n , so that s is higher than L + A/2 during non negligible time intervals 
(dashed lines represent s = Dsi n ) in) during such a time intervals L — mm k (S^.(q k )) 
(or inhij(Sj(yj))) is increasing non negligibly towards s, so that L cannot both 
converge towards L and stay lower than L during the whole time interval: there is 
a contradiction. 



Substep 3b. 2: s is higher than L+X/2 during non negligible time intervals. 

Because of the boundedness of s 

s < Ds in 

every time s is higher than L + X, it has been higher than L + X/2 during a non 
negligible time interval of minimal duration A(X) — 9 p s . On Figure |#| the case 
s = Dsi n is represented by dashed lines. 

Substep 3b. 3: L = rnin^jSf (cjfc)) (or mm.j(Sj(yj))) is increasing non negligi- 
bly towards s, so that L cannot both converge towards L and stay lower 
than L during the whole time interval: there is a contradiction. 

Like in previous proofs, we are interested in what happens during the [t s — A(X),t s ] 
time-interval, with t s — A(X) > t L (e) (for some e < X). Since, during this time- 
interval, s(t) > L + X/2 and L < L, we know that there exists a k such that 
L(t s ) = S£(q k (t s )) <L, or a j such that L(t s ) = S^(y 3 (t s )) < L. 

For both this step (3b. 3) we choose to first only present arguments for the case 
L(t s ) = mmk(S^,(qk)); almost similar arguments for the case L(t s ) = minj(S'J (yj)) 
will then be briefly presented. 

• if L(t s ) = mmk(Sj:(qk)), then during the whole considered time-interval, as 
S^{qk) was increasing, we know that 

L-e<L< S z k {q k ) < S z k {q k {t s )) < L (27) 

so that Q k (L — e) < q k (t) < Q k (L). For the k species, the dynamics of q k can 
then be lower bounded: 

q k >Pk(L + X/2)-f k (Q k (L)) 
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and then 

q k > p k {L + X/2) - p k (L) = G k (X) 

positive, so that the increase of q k during the [t s — A(X),t s ] time-interval is 
also lower bounded: 

q k (t s ) - q k (t s - A(X)) > G k (X)A(X) = H k (X) 

Since Q k = S k is locally Lipschitz with constant K (because f k > 0), we 
have 

q k (t s ) - qu(t s - A(X)) = Q k (S z k (q k (t s ))) - Q k (S z k (q k (t s - A(X)) j) 
<K[S* k (q k (ts))-S* k (q k (ts-A(X)))] 

so that the corresponding increase of S k (q k ) is lower bounded with 

S* k (q k (t*)) - S* k (q k (t* - A(X))) > ±H k (X) 

and then 

S* k (q k (t s -A(X)))<L-±H k (X) 
which implies the same higher bound for L: 

L(t s ~A(X))<L-^-H k (X) 



By choosing e < j^H k (X), this inequality is contradictory with (27) so that 
Case 2 is not possible 

• if L(t s ) — mh\j(S V j{yj)), then the same arguments can be developped for the j 
species, with a lower bound Gj(X) on the yj dynamics: 

Gi (A) = ft- (L + A/2, Y, (L)) - fa (L, Y 3 (L)) > 

and then an increase of variable yj at least equal to Hj(X) = Gj(X)A(X) 
followed by a non negligible increase of L 

L(t s — A{X)) < L — j^Hj(X) 

because Yj is locally Lipschitz. Finally a contradiction also occurs when e < 
^H 3 {X): 

L(t s -A(X)) <L-e 

C Computation of system £ Jacobian Matrix and 
eigenvalues for all the equilibria 

Computation of the Jacobian Matrix of system S, with s = s;„ — x i~~ Sfci Vi~ 

l^ k= l QkZk- 

/ j xx j x y j xz j xq \ 

jyx jyy jyz jyq 

jzx jzy jzz J zc i 
\ J qX J qy J qZ J qq J 
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where 



Jff = ai (s) - D - and VI ^ i, J™ = -^ Xi 




jyx _ OPj 

J]I=PAs,y j )-D+^y j -^y j and ^j,Jjf = -f%- 

jyz 3f5j 

J jk ~ ~l)s~yj1 k 




and 




Jul = lk{q k ) - D and Vl^k,J z k f=0 
Kl = ltzk and Vl^k,J z k ? = 



and 

jqx _ dp k 

J ki — ds 

jIV _ dpk 

J kj — ds 



Fortunately for eigenvalue computations, at equilibria the null biomasses will 
simplify the matrix: 



• when Xi = 0, then the whole i th line gives eigenvalue cti(s) — D (denoted 
"xi- eigenvalue") and can be deleted, as well as the i th column; 



• whenyj = then the whole N x +j th line gives eigenvalue f3j(s,yj)—D (denoted 
"yj -eigenvalue") and can be deleted, as well as the N x + j th corresponding 



• when Zk = then the whole N x + N v + k th line gives eigenvalue jk(q k ) — D 
(denoted "z^- eigenvalue") and can be deleted, as well as the N x + N y + k th 
column; in a second step, the whole N x + N y + N z + k th column can also be 
deleted and gives eigenvalue (denoted "q^- eigenvalue") , as well as the 

N x + N y + N z + k th line. 

C.l Complete washout equilibrium 

With this in hand, we see that for equilibrium Eq (x^ = yj = = 0) the Jacobian 
matrix is triangular, so that the eigenvalues lay on the diagonal. They are: 



column; 



• Oii{s in ) - D 



• Pj(8 in ,0)-D 

• lk{Qk{s m )) - D 

• — 5^" (negatives) 
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We denote n x , n y , n z the number of M-, C- and Q- species verifying the inequalities 
of Hypothesis [5[ and thus having the possibility to be at equilibrium with a positive 
biomass, under controls D and Sj n . Each of these species has a positive correspond- 
ing eigenvalue on this equilibrium, so that equilibirum Eq has n x +n y + n z positive 
eigenvalues, and N x — n x + N y — n y + 2N Z — n z negative eigenvalues. 



C.2 M-only equilibria 

For equilibrium E'f we get all the previously cited x-,y-, z-and q-eigenvalues: 

• ai(sf*) — D whose signs are the same as sign(sf* — sf*) 

• j3j(sf*, 0) — D which are positive if the j species is sf* -compliant, or negative 
else; 

• lk{Qk{sf*)) — D whose signs are the same as sign(sf* — s|*) 

• ~Qq k k which are all negative 

and the remaining eigenvalue corresponds to the positive Xi-only dynamics: 

X% xf) D)xi 

which yields the eigenvalue —^§^x* for free bacteria species i. Each free species 
with a substrate subsistence concentration sf* or s^* lower than sf* gives a positive 
eigenvalue. Among all the E'f equilibria, only Ef is stable if and only if s* = sf* < 
sf* , and if all the C- species are not sf* -compliant. 



C.3 Q-only equilibria 

For Equilibrium Ef, we get all the 

• x-eigenvalues whose signs are the sign of sign(sf,* — sf*); 

• y- eigenvalues: as previously, y- eigenvalues are positive if the corresponding 
C-species is sf? -compliant and negative else; 

• z\- eigenvalues whose signs are the sign of sign(sf* — sf*); 

• qi- eigenvalues for all I ^ k (negative); 

and the remaining eigenvalues correspond to the positive (zk, qk)-only dynamics: 

h = (lk(qk)-D)z k 
qk = Pk{sin - qkZk) - fk(qk) 

and we obtain the following resulting matrix: 

( o° o^ Zk 9f ) 

\ ds y fc ds 4k dq k J 

which has negative trace and positive determinant, so that its two eigenvalues are 
real negative. Just like before, each free species with a substrate subsistence con- 
centration sf* or sf* lower than sf* gives a positive eigenvalue. Among all the Ef. 
equilibria, only Ef is stable if and only if s* = sf* < sf* , and if all the C-species 
are not sf* -compliant. 
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C.4 C-only equilibria 

Now let us consider the E G equilibria for which all j G G (where G represents a sub- 
set of {1, . . . ,N y }) C-species coexist in the chemostat under substrate concentration 
s V q , while all the free species are washed out. s V q is defined by s G + X}jeG ^j( s G ) = 
s in . Note that some of the G species can have a null biomass on these equilibria, as 
Yj(sq) might be null for some j G G. 
This gives all the 



x-eigenvalues whose sign are the same as the signs of sf* 
z-eigenvalues whose sign are the same as the signs of si* 



„v*. 
„v*. 

S G > 



• q-eigenvalues (negative). 

All the yj species who are not included in G give negative eigenvalues if they are 
not Sq -compliant, and positive eigenvalues else; their eigenvalues cannot be null 
because of technical hypothesisVA All the yj species who are included in G but have 
a null biomass Yj(s v G ) on the Eq equilibrium give negative eigenvalues. Now let us 
study the remaining matrix J G V which is composed of all the j G G lines of J vv , for 



which Yj(s v G 



> 0, and thus (3j(s G 



Y j( s 



V g)) 



D: 



fit 



which yields the Jacobian matrix: 

( -a x -bi ... -ai 



jyy _ 



-aj - bj 



-a n - b n J 



with aj = ^-Yj(s V Q) > and bj = -^Yj{s v G ) > 0. 

Let us show that this matrix has only real negative eigenvalues, by using the 
definition of an eigenvalue A = (A + Bi), where A € M is the real part and Bel 
the imaginary part: 





/ 2/i > ^ 




( yl \ 


ryy 
• J G ' 


\Vn J 


= (A + Bi) 


\Vn J 



(28) 



We obtain n equations: 



id thus 



-bjUj - aj ^2vi = (A + Bi)yj 
i 

(A + Bi + bj) Vj = -aj Vi 



(29) 



RR n° 8038 



40 



Masci, Grognard, Benoit & Bernard 



If we have A + Bi + bj =0 for some j , then B — and A = —bj < so that we 
have a negative eigenvalue. 
Else, isolating yj yields 



bj + A + Bi 



Summing over j , we obtain 



3 



- a 3 Ej vi 

i n + A + Bi 



Now i/Ej Uj ~ 0? since some yj must be different of 0, (29) yields, for that j, that 
A + Bi + bj =0 so that again B = and A = —bj < 0. 
Else, simplifying the sums of yi andyj, this yields 

Ej V bj+A+Bi 

-a.j (bj+A—Bi) 

,B) 



Ej \ (b 3 +A) 2 + B 



Since the left-hand- side is real, the imaginary part of the right-hand side must be 
zero, which imposes B = 0. For thr right-hand-side to be positive, at least one of 
the bj + A must be negative, which translates into mirij (bj + A) < and 



A < — min 6j < 

j 

We conclude from this that all eigenvalues of this matrix are real negative. 

Finally, an E V G equilibrium is stable if and only if all the C-species not contained 
in G are not s V q -compliant (this is equivalent to saying that s V q = s y * , with s y * = 

S {l,...,Ny})> ^ it 3 * ^ SV *- 



C.5 M-coexistive equilibria 

In this section we consider equilibria e[ X q^ where free bacteria species Xi coexists 
with the C-species in G, a subset of {1, . . . , N y }, under substrate concentration sf*. 
We obtain here all the 

• x\- eigenvalues (I ^ i) whose signs are the signs of sf* — sf*; 

• z-eigenvalues whose sign is the sign of sf* — sf,* ; 

• q-eigenvalues (negative); 

yj -eigenvalues with j not in G are positive if yj is sf* -compliant and negative else; 
yj-eigenvalues with j in G but have a null biomass Yj(sf*) give negative eigenvalues. 
For the remaining C-species, and species Xi, we obtain the following system: 

( ±i = (ai(s in - Xi-J^iVl) ~ D)xi 

\ in = Wj ( s in - xi -Y.iVuvs) - D ) V] 
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and the Jacobian matrix: 



with ao 

j e {i, • 



-a 
-ai 



—a 
-oi - 61 



-«o 



*j - b j 



\ -<h> 



-a 
-«1 



6„ / 



^< > 0, a 3 = ^{sf ) > and b, = -g^OT) > <V 
, n} ). This matrix has exactly the same form has the one considered on 
Appendix \C4\ The only difference being that the there is no "bo" in the first element 
of the matrix. Defining a bo — 0, we can then conclude that all eigenvalues are real 
and negative because, following the development of Appendix \ C.j\ we obtain 

A < - min b, = 



Finally, only equilibrium E^ l x 



{l,...,N y } 



can be stable if and only if s* = sf* 



C.6 Q-coexistive equilibria 

In this section we consider equilibria E^q^ where phytoplankton species Zk coexists 
with the attached species in G, a subset o/{l, . . . , A^}, under substrate concentration 

b k ■ 

We obtain here all the 



• x-eigenvalues whose signs are the signs o/s|* 

• zi- eigenvalues (I 7^ j) whose sign are the signs o/s|* — sf*; 

• q-eigenvalues (negative); 

yj -eigenvalues with j not in G are positive if yj is sf* -compliant and negative else; 
yj -eigenvalues with j in G but have a null biomass Yj(st*) give negative eigenvalues. 
For the remaining C-species, and species z k , we obtain the following model 

= (Pj(sin -J2iVl~ IkZkiVj) - D)yj 
= (ik(qk) - D)z k 

= Pk{sin - YsiVi - ikZk) - fk{qk) 

and, swapping the last two equations and using fk{qk) — lk{qk)ok, we get the Jaco- 
bian matrix: 



( -ai - 61 



-ax 



-at 







>>n+l 



\ 



-a\Zk 



-CljZ k 



-aiq k \ 



-a 3 q k 



0"n Un 



-a n+ iZk — b n+ iqu — 7 — a n +iqk 

bn+lZk / 
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with a 3 = ^Y 3 {s z *) > and b } = -fff^'Of) > f or 3 £ {V" ,n}, with 
a n+i = ^§f~ an d b n +i = f^-- By using the definition of eigenvalue A = (A + Bi) (see 
(28)) we follow a similar path to that of Appendix G4- we show that the eigenvalues 
are real and negative. 



Finally, only equilibrium E^ z ^ N j can be stable if and only if s* = sf 

Remark 6 The same work can be done for the whole system where the eigen- 
values are the same, plus the —D eigenvalue which arises from mass balance dy- 
namics 13). 
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